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Abstract 

We report results of molecular dynamics simulations of liquid water at the temper- 
ature T = 277 K for a range of high pressure. One aim of the study was to test the 
model Amoeba potential for description of equilibrium structural properties and 
dynamical processes in liquid water. The comparison our numerical results with the 
Amoeba and TIP5P potentials, our results of ab initio molecular dynamics simu- 
lations and the experimental data reveals that the Amoeba potential reproduces 
correctly structural properties of the liquid water. Other aim of our work was re- 
lated with investigation of the pressure induced structural transformations and their 
influence on the microscopic collective dynamics. We have found that the structural 
anomaly at the pressure pc ~ 2000 Atm is related with the changes of the local, 
short-range order in liquid water within first two coordination shells. This anomaly 
specifies mainly by deformation of the hydrogen-bond network. We also discuss in 
detail the anomalous behavior of sound propagation in liquid water at high pressures 
and compare numerical results with the experimental data. 

Key words: Amoeba water model, molecular dynamics, structural anomaly, 
dynamic structure factor, hydrogen-bond network 
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1 Introduction 



Water is one of tlie most widespread liquids in wildlife and has a great num- 
ber chemical and technological applications [1]. Although a particular water 
molecule has a simple chemical structure, a water system is considered as a 
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complex fluid because of its anomalous behavior in tliermodynamical, struc- 
tural and transport properties [2]. It dilates with solidification, and the density 
has a maximum at the pressure 1.0 Atm and the temperature 277 K [3]. Ad- 
ditionally, there is a minimum in the isothermal compressibility at 319 K and 
a clear minimum in the isobaric heat capacity at 308 K. These anomalies are 
linked with the microscopic structure of liquid water [4], which can be re- 
garded as a transient gel - a highly associated liquid with strongly directional 
hydrogen bonds [5,6]. 

An active study of the different characteristics and properties of water system 
is provided by the development of the numerical methods of the computer 
simulations. So, the first study of the water "structure" was performed by 
Barker and Watts within Monte-Carlo method [7]. They have computed the 
energy, heat capacity and radial distribution function for water system at the 
temperature 298 K and have compared the obtained results with the experi- 
mental data. The considered system was very small, 64 molecules interacted 
via the Rowlinson potential [8]. In the next study performed by Rahman and 
Stillinger [9] the water system was simulated through 216 rigid molecules in- 
teracted by means of the pair-additive effective potential. The authors have 
investigated in detail the structural, transport and dynamic properties of liq- 
uid water at the temperature 307.3 K. The molecular dynamics simulations in 
the low temperature thermodynamic phase range were performed by Ref. [10]. 
Here, the "second critical point" at the temperature T = 223 K and pressure 
p — 1000 Atm was found. Below this second critical point, the liquid phase 
separates into two distinct phases - a low-density liquid and a high-density liq- 
uid. It is necessary to note that to study the water system a great number of 
the model potentials for intramolecular and intermolecular interactions were 
suggested (sec, for example, review [11]). Here, the first model of the liquid 
water was proposed in 1933 by Bernal and Fowler [12]: an ice-like disordered 
tetrahedral structure arising from the electrostatic interactions between close 
neighbors. 

The special role in the understanding of the microdynamic features of the wa- 
ter have the methods of ab-initio simulations based on the density functional 
theory. So, one of the first works, where the Car-Parrinello method was ap- 
plied to study dynamic properties in water, was done by Laasonen et al. [13]. 
They also found that the local density approximation gives realistic results 
for the intramolecular properties without gradient corrections, but it fails to 
give a correct description of the intermolecular interactions and the binding 
properties. Later, the numerous ab-initio simulations for water system have 
been done [14,15,16,17,18]. 

Recently, a large amount of the works were intended for the investigation of 
water structural, transport, mechanical properties at the different externally 
apphed conditions [19,20,21,22,23]. For example, the processes of the disso- 
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ciation of the water molecules at the pressures 14.5 GPa and 26.8 GPa were 
investigated in Ref. [24] by means of ab-initio molecular dynamics simulations. 
Molecular dynamics simulations of water at the negative pressures were car- 
ried out with two model potentials of interparticle interaction: the extended 
simple point charge (SPC/E) and the Mahoney-Jorgensen transferable inter- 
molecular potential with five points (TIP5P). The equilibrium phase diagram 
was constructed in Ref. [25] for a wide range of temperature and pressure by 
means of Monte-Carlo simulations within two model potentials, TIP4P and 
SPC/E. Then, the correlation between structural and dynamical anomalies in 
supercooled water was studied [26] on the basis of the SPC/E potential. 



Some of works were pointed mainly towards the study of structural and 
dynamical properties related with the anomalous features of liquid water 
[27,28,29,30,31,32,33,34,35]. The concepts of low-density water (LDW) and 
high-density water (HDW) were suggested in Ref. [27] on the basis of the 
experimental data of neutron diffraction, which detected clear the structural 
transformations in liquid water [28]. Further, the structural transition from 
LDW to HDW as well as the anomalies related with these structural trans- 
formations were extensively considered [29,30,31,32,33,34]. An comprehensive 
overview of the properties of supercooled and glassy water, the transition be- 
tween low-density amorphous (LDA) ice and high-density amorphous (HDA) 
ice were presented in Ref. [35]. Nevertheless, in spite of the numerous studies, 
the question about the impact of structural transition from LDW to HDW on 
the dynamical properties and anomalies of water is still open [3]. 



In this work we perform the numerical study of the structural transformations 
in water, which are appeared at the fixed temperature T = 277 K with the 
change of the externally applied pressure p = 1.0 ^ 10 000 Atm. In addition, 
the influence of the structural rearrangements on the dynamical properties in 
liquid water is also considered. Other aim of the work is related with the test 
of the Amoeba potential [36,37,38] to describe the structural and dynamical 
features of water. 



The paper is organized in the following way. In the next section, we describe 
the details of the molecular dynamics simulations of liquid water on the basis 
of the Amoeba model potential. A detailed analysis of the pressure dependence 
of the static properties and the structure anomalies are then reported. Then, 
the features in the dynamic and dispersion of sound velocity are reported and 
compared with the experimental data on inelastic X-rays scattering. The role 
of the hydrogen-bond network in the structural transformation is discussed. 
And, finally, we conclude with a short summary. 
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2 Computational Details 



Equilibrium classical molecular dynamics (MD) simulations of water at the 
constant temperature T — 277 K and the different pressures were carried 
out with the TINKER molecular modeling package [36]. We use the Amoeba 
(Atomic Multipole Optimized Energetics for Biomolecular Applications) water 
model, which was recently suggested by Ren and Ponder [37,38]. Our compu- 
tations were performed for 4000 water molecules interacted within a cubic 
box, where the periodic boundary conditions were imposed in all directions. 
The Ewald summation was used to handle the electrostatic interactions as 
well as an atom-based switching window of the size 12 A was applied to cut 
off the van der Waals interactions. The equations of motion were integrated 
via a modified Beeman algorithm [39] with the time step Ar = 1.0 fs. The 
MD calculations were performed in the isothermal-isobaric ensemble at the 
temperature T = 277 K and the pressure p = 1.0, 1000, 2000, 3750 5000, 7797 
and 10000 Atm, where we used the Berendsen thermostat and barostat [40] 
to enforce the constant temperature and pressure. 



3 Results 

3.1 Partial Radial Distribution Functions 

The radial distribution function (RDF) provides directly the information about 
the average packing the molecules in water and can be extracted from neu- 
tron [41] and X-ray scattering data [42]. On the other hand, this term in an 
extended form can be calculated on the basis of MD simulation data [43], 
namely 



where ga,p{f^) is the probability to find an atom in the range [r, r -|- Ar] and L 
is the length of the box edge. The indexes a and /3 define the type of the atom, 
a, /9 G {O, whereas and iVg are the particle numbers of the type a and 
/3, respectively. The quantity nji3{r) specifies the number of ^-atoms located 
at a distance r from a jth atom within a spherical layer of the thickness Ar. 

In Fig. 1, the partial RDF's of atoms 5'a,/3(r) {a, (3 G {0,H}) obtained from 
our simulations with the Amoeba potential at the temperature T = 277 K 
and the pressure p — 1.0 Atm are compared with the molecular dynam- 
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Fig. 1. (color online) Partial radial distribution functions of atoms ga^pi^) 
{ct,/3 £ {0,H}) at the temperature T = 277 K and the pressure p = 1.0 Atm. 
Solid line represents the results of classical molecular dynamics simulation with 
the Amoeba potential; dotted line show the results of molecular dynamics simula- 
tions with the TIP5P potential; dashed line represents the results of Car-Parrinello 
molecular dynamics simulation for the system with 64 molecules [47] ; circles are the 
experimental data on neutron diffraction [48]. 



ics simulations with the TIP5P potential, results of ab-initio Car-Parrinello 
molecular dynamics [47] and experimental data on neutron diffraction [48]. 
Although molecular dynamics simulations data with the Amoeba potential 
overestimate insignificantly the high of the first peak in partial RDF's in com- 
parison with experimental data and results of ab-initio molecular dynamics 
simulations, they reproduce correctly the structure of liquid water at the tem- 
perature T = 277 K and the pressure p = 1.0 Atm. Moreover, the comparison 
reveals that the model TIP5P has a better agreement with experimental data 
and results of ab-initio molecular dynamics simulations. 

The pressure dependence of the partial RDF for water is presented in Fig. 2. 
As one can see from the figure, the partial contributions of RDF for the pairs 
H-H and 0-H remain unchanged practically at the increase of the pressure. 
The remarkable changes are observed in the pressure dependence of RDF for 
0-0 pairs and indicate clear on the transition pressure ptr ~ 2000 Atm. It 
is interesting that the atomic rearrangements corresponding to the structural 
transformation are restricted to the second coordination sphere, where the 
shift of the second peak in goo{r) to the range of the high r occurs. 
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Fig. 2. (color online) Partial radial distribution functions of atoms Qa^jiiT) 
(a, /3 G {O, if}) at the temperature T = 277 K and the different pressures. 
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3.2 Partial Static Structure Factors 



Other term, which is sensitive to the structural transformations, is the static 
structure factor, S{k). This term is measured experimentally and is defined 
as [43] 

TV , N 



1 / ^ 1 ^ 

\j=l 1=1 



e^^^ ' . (2) 



Here, k is the wave-vector and k = |k|; rj(t) defines the radius- vector of an jth 
atom. The partial static structure factor Sa,i3{k), as an extension of S{k), is 
directly related with the Fourier transform of the radial distribution function 
[39]: 



S, 



a,fi{k) = l + Ann j 





9a A'^) - 1 



sin(/cr) 
kr 



dr.. 



(3) 



where n — N/V is the numerical density. 

The partial static structure factor Sa^pik) at various pressures for different 
atom pairs, H-H, 0-H and 0-0, is presented in Fig. 3. As can be seen, the 
features in the first peak of Sa,p{k) become more pronounced with the in- 
crease of pressure that indicates on the rise of the short-range order. Further, 
the split of the first peak in Soo{k) disappears with the increase of pressure 
that is related with the local molecular ordering in the vicinity of the first 
coordination shell. 



3.3 Wendt- Abraham Parameter, Pair- Correlation Entropy and Translational 
Order Parameter 



An additional information about the local structural features of the system 
can be extracted with the help of the extended Wendt- Abraham parameter 

„rnin( \ 

WA ^ 9a,p 

9a,p V) 



where g^1f{r) and g^'^^{r) represent the main maximum and the first mini- 
mum in a partial radial distribution function, respectively [49]. The thermody- 
namic excess entropy of a fluid is deflned as the difference in entropy between 



7 




2 000 



4 000 6 000 8 000 10 000 
Pressure (Atm) 



Fig. 4. (color online) Pressure dependence of the partial Wendt- Abraham parameter 
r^^, the partial pair-correlation entropy 82'^ and the partial translational order 
parameter Pa^p, where a,l3 G {0,H}. 

the fluid and the corresponding ideal gas under identical temperature and 
density conditions. The total entropy of a classical fluid can be written as 



N 



s = s, 



id 



71=2 



(5) 



where Sid is the entropy of the ideal gas reference state, Sn is the entropy 
contribution due to n-particle spatial correlations. Then, the excess entropy 
is deflned as 

Sex = S — Sid- (6) 
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The main contribution of S^x for classical fluids is a pair-correlation entropy 
5*2 (for example, in case of a monoatomic liquids this two-particle contribution 
to the excess entropy is 85 -j- 95% over a fairly wide range of densities), which 
can be given in the extended form as 



oc 

S^'^ = -27rn j [g^,p{T) ln(^„,^(r)) - [g^,p{r) - Ijjr^dr, (7) 



where the term 5*2 is measured in units oiks- 

The translational order parameter p is a simple and convenient measure to 
test the radial ordering in the studied system [50]. The separation of p into 
partial contributions can be written as the following 



Poc,p = — / \9ocA^) - ^\dr, (8) 



where is the distance at which the structural features of the system cease 
to be considered. The value of corresponds to a half of the edge of the 
simulation box. It is clear that an increase of Pa^p indicates directly on the 
growth of structural radial ordering, whereas the increase of r^J- and 82'^ 
manifests the processes of structural disordering in the system. 

Figure 4 illustrates the pressure dependence of the partial Wcndt-Abraham 
parameter r^^, the partial pair-correlation entropy 82'^ and the partial trans- 
lational order parameter p^^^ at the fixed temperature T — Til K. The de- 
crease of indicates on the increase of the two-particle chain clusters in 
the system, which appear, e.g., due to the system densification and precedes 
the transition of the system into a solid (crystalline or amorphous) phase. 
As can be seen from the figure, all parameters discover the change in the 
pressure dependence at Pc ~ 2000 Atm. It is remarkable that the extended 
Wendt- Abraham parameter is more sensitive to detect these structural 
transformations than other parameters, 82'^ and p^,^- This can be consid- 
ered directly as an evidence that the structural transformations detected at 
Pc ~ 2000 Atm are related with the changes of the local, short-range order, 
since the parameter probes the structural properties inside of the first 

coordination shell whereas the quantities 82'^ and p^,^ characterize features 
at more extended spatial scales. 
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Fig. 5. (color online) Main: Tetrahedral order parameter versus pressure at the 
temperature T = 277 K. Inset: Instantaneous configurations of water molecules 
at the two values of the tetrahedral order parameter: = 0.01 (the four bonds 
are practically aligned) and Qk = 0.99 (molecular system is close to be a perfect 
tetrahedron) . 

3.4 Tetrahedral Order Parameter 



The tetrahedral order parameter can be used to quantify the tetrahedrality 
within the first shell [28]. This parameter is defined as 



q 3 4 



COS 6ikj + 



(9) 



The term Oikj is the angle formed between the central molecule k and its 
neighbors i and j. The average value 

1 ^ 

{Q) = i;7EQ^ (10) 

k=l 



quantifies directly the orientational order within the first shell. For the perfect 
tetrahedral order, one has (Q) = 1; if the bonds are arranged in a manner with 
the absence of tetrahedral geometry, one has (Q) = 0. So, with the deviation 
from the perfect tetrahedron the values of (Q) demonstrate decease. 
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The pressure dependence of the tetrahedral order parameter (Q) for hquid wa- 
ter at the temperature T = 277 K is presented in Fig. 5. As can be seen, with 
the increase of the apphed pressure the tetrahedral order parameter demon- 
strates a nonhnear decrease. Such a behavior of the parameter indicates on 
the decrease of the amount of tetrahedral structures in liquid water with the 
isothermal increase of the pressure p. that is in agreement with the results 
of Ref. [27]. It is necessary to note that although the behavior of the tetra- 
hedral order parameter (Q) changes with the pressure p, this parameter (Q) 
demonstrates no a sharp transition from low-density phase to high-density 
one. 



3.5 Dynamic Structure Factor and Dispersion of Sound Velocity 



According to a definition (see Ref. [43]), the dynamic structure factor S{k,u) 
is related with the Fourier transform of the intermediate scattering function 
F{k,t) 

S{k, uj) = - f F{k, t) cos{ujt)dt. (11) 

TT J 




Here 



Fik,t) 




ik-Viio) 



3 



-jk.r,(t) 



{5p,{Q)5pl{t)), 



[12) 



where 6pk{t) is the fc-component of the fluctuation of the microscopic number 
density pk{t). Nevertheless, because of the long-time tails of F{k,t) in water 
(see Ref. [51]), the finding of the dynamic structure factor spectra S{k,u) 
by means of Eqs. (11) and (12) yields the significant inaccuracies. To avoid 
this problem we use the Wiener- Khinchin (statistical-autocorrelation) theorem 
which allows one to rewrite Eq. (11) in the next identical form 



S{k,uj) = — 
tu 



6pk{t)e"^'dt 



(13) 



where tM is the observation time for the variable Spk{t). Eq. (13) gives a 
possibility to find S{k,uj) from the dynamical variable Spk{t), directly. Note 

that an advantage of Eq. (13) over Eq. (11) is related with the saving of 
computational time, that is very important at numerical estimations. Then, 
the dynamic structure factor S{k, cu) can be evaluated at these /c-values, which 
are allowed by the system size: k — 27r/L(n, m, I), where n, m, I are integers. 
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Fig. 6. (color online) Dynamic structure factor of liquid water at the temperature 
T = 277 K and the pressure p = 3750 Atm for different values of the wave number 
k: the noisy (gray) hues represent the results obtained from molecular dynamics 
simulation data and Eq. (13), thick (black) lines are results of the smoothing pro- 
cedure. 

The size of our system allows one to take the smallest value oik ^l.?>% nm~^. 
To delete the noises from dynamic structure factor spectra S{k^ u) we apply 
the smoothing procedure based on the computation of the arithmetic mean 
for the local window with the length 1.0 ps~^. 

The frequency spectra of the dynamic structure factor S{k,uj) of water for 
different values of the wave number k = 2.69 ^ 8.63 nm~^ at the temperature 
T = 277 K and the pressure p = 3750 Atm (that is higher than p^) are 
presented in Fig. 6. Here, the high frequency collective excitations responsible 
for sound propagation in the system are clearly detected for all the considered 
values of k. To obtain the precise values of the sound velocity u!c{k) which 
is related with the high-frequency peak-position in S{k,uj), we consider the 
longitudinal current spectra CL{k,uj) [53]. This term can be found from the 
dynamic structure factor S{k, cu) by the way: 



CL{k,uj) = ^S{k,u). 



(14) 



On Fig. 7, the dispersion of sound velocity ujc{k) in water at the temperatures 
T — 277 K and T = 297 K and the pressures p = 3750 Atm and p ~ 
7797 Atm are presented. Results of molecular dynamics simulations with the 
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Fig. 7. (color online) Main: Dispersion of sound velocity ujc{k): circles represent 
the experimental data [54]; triangles are results of molecular dynamics simulations 
data with Amoeba potential; squares represent results of computer simulations with 
TIP5P potential; solid line indicates the best linear fit to experimental data, and 
the slope corresponds to the sound velocity "i? = 3320 m/s {'d = 3540 m/s) at the 
temperature T = 277 K (T = 297 K). Inset: Partial dispersion of sound velocity 

Amoeba model potential reproduce correctly the experimental inelastic X-ray 
scattering data [54]. Moreover, with the increase of k the simulations results 
overestimate slightly the experimental data for sound velocity of Ref . [54] . The 
comparison of the results with the two potentials reveals that data obtained 
on the basis TIP5P are in the range of inaccuracies of the dispersion curve 
obtained on the basis Amoeba potential and have a better agreement with 
experimental data. Nevertheless, simulation data with TIP5P demonstrate 
also overestimated values in comparison with the experimental outcomes of 
inelastic X-ray scattering. 

The partial dispersion of sound velocity uj"'^{k) for O and H atoms obtained 
with the Amoeba potential is given in inset of Fig. 7. As can be seen, both 
lines are identical. This is direct evidence that the main influence on the 
features of sound propagation in water appear due to oxygen, which is the 
heavy component of a water molecule. Note that this is agreed with results of 
Ref. [55] (p. 7235). 

To estimate quantitatively the pressure dependence on the properties of the 
collective dynamics in water we show in Fig. 8 the density and pressure depen- 
dence of the sound velocity where results of our molecular dynamics simu- 
lations data are compared with experimental data [54]. First, as can be seen 
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Fig. 8. (color online) Main: Density dependence of the sound velocity i?: the circles 
represent the experimental IXS data [54], the triangles are results of our molecular 
dynamics simulations with Amoeba water model. Solid line is the spline-interpola- 
tion. Inset: Pressure dependence of the sound velocity obtained from simulations. 



from figure, results of the molecular dynamics simulations reproduce correctly 
the experimental data for sound velocity [54]. An insignificant understatement 
of i) appears from simulations in comparison with experimental data at high 
densities. Nevertheless, both our numerical results and experimental inelas- 
tic X-ray scattering data reveal clearly a changes in dynamics at the density 
Pc = 1.114-f- 1.12 g/cm^. Note that such density of the investigated water sys- 
tem corresponds to the pressure pc ~ 2000 Atm (see inset of Fig. 8), at which 
the structural transformations were detected by the order characteristics. 

Thus, the structural and dynamic anomalies observed in water at the pressure 
Pc ~ 2000 Atm are provided by the local structural molecular rearrangement 
inside of the first two coordination shells. It is necessary to note that such a 
treatment is also suggested by results of Ref. [32], where anomalies in water 
were studied on the basis of the five-site transferable interaction potential 
(TIP5P). It was found in Ref. [32] that the anomalous decrease of orientational 
order upon compression occurs in the first and second coordination shells, but 
the anomalous decrease of translational order upon compression occurs mainly 
in the second shell. 
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Fig. 9. (color online) Pressure dependence of the hydrogen bond angle. 
3. 6 Hydrogen-Bond Network 

A detailed description of the hydrogen-bond network in liquid water is the key 
to understanding its unusual properties. The studies of the hydrogen-bonded 
network structure in water have mainly relied on neutron and X-ray diffrac- 
tion, infrared spectroscopies, and computer simulations techniques [56]. Ex- 
perimental data from noncrystalline materials provide the radial distribution 
functions and static structure factors [48] that do not provide angular corre- 
lations needed to assign uniquely the local geometries in water [57]. In order 
to understand the effect of the structure on the dynamics we carried out a de- 
tailed analysis of the pressure dependence of the hydrogen bond angle [58,59]. 
In the same lines with the work [60], we use a geometric criterion. According 
to this criterion two water molecules are taken to be hydrogen-bonded if their 
inter-oxygen distance is less than 3.5 A and simultaneously the hydrogen- 
oxygen distance is less than 2.45 A and the oxygen-oxygen-hydrogen angle 
/3 is less than 30°. We note that the critical distances of 3.5 A and 2.45 A 
are essentially the positions of the first minimum in the oxygen-oxygen and 
oxygen-hydrogen radial distribution functions, respectively. All the terms of 
this algorithm are presented in the bottom inset of Fig. 9. The behavior of 
the hydrogen bond angle {3 at the temperature T = 277 K is illustrated in 
Fig. 9. This quantity is presented as an averaged over all possible hydrogen- 
bond environments. As can be seen from the figure, the distinct changes in 
the behavior of the angle (3 at critical pressure p = 2000 Atm is observed. 
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In the inset of Fig. 9, wc show the configurations of water molecules at the 
applied pressure p = 1.0 Atm (on the left of Fig. 9) and p = 10000 Atm (on 
the right of Fig. 9) in the neighborhood of four closest molecules, which are 
connected by hydrogen bonds. Thereby, we confirm quantitatively the idea 
that the structural and dynamical anomalies of water are related with the lo- 
cal structural rearrangement of molecules within first two coordination shells 
[32], which specifies by deformation of the hydrogen-bond network. 



4 Conclusions 



Molecular dynamics simulations of liquid water at the temperature T = 277 K 
and the range of pressure p = 1.0 -j- 10 000 Atm are performed on the basis 
of model Amoeba potential [36]. It is found that this potential reproduces 
correctly equilibrium structural and dynamical properties of liquid water, al- 
though TIP5P potential yields a better agreement with experimental data. 

The pressure dependence of the structural characteristics and order parame- 
ters detects clearly the appearance of the structural anomaly [32] in the water 
system at the temperature T — 277 K and the pressure p ^ 2000 Atm. The 
detailed analysis of the extended Wendt- Abraham parameter, pair-correlation 
entropy and translational order parameter indicates that the observed struc- 
tural transformations occur in the hydrogen/oxygen atoms in the range of the 
first two coordination shells. Thus, these structural transformations are cor- 
related with the deformation of the hydrogen-bond network [59] . Moreover, it 
is found that with the increase of the applied pressure, a fraction of tetrahc- 
dral structures in the system decreases. The tetrahedral order parameter (Q) 
demonstrates no a sharp transition from low-density state to high-density one. 

We have found that the structural anomaly has an impact on the dynamic 
properties of the liquid water. The computed spectra of the dynamic struc- 
ture factor S{k, u) at high pressures reveal all features (including the high- 
frequency acoustic mode) typical for liquid systems [4,61,62]. The found sound 
dispersions are in a good agreement with the experimental inelastic X-ray scat- 
tering data [54]. It is found the main features of sound propagation in liquid 
water are related with the oxygen atom. 

The pressure dependence of the sound velocity indicates on the anomalous 
dynamics in water at the same value of the pressure as for the structural 
anomaly, p ~ 2000 Atm (the temperature T = 277 K), that corresponds to 
the density Pc ~ 1-12 g/cm^. 
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